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We model recruitment in adaptive social networks in the presence of birth and death processes. 
Recruitment is characterized by nodes changing their status to that of the recruiting class as a 
result of contact with recruiting nodes. Only a susceptible subset of nodes can be recruited. The 
recruiting individuals may adapt their connections in order to improve recruitment capabilities, thus 
changing the network structure adaptively. We derive a mean field theory to predict the dependence 
of the growth threshold of the recruiting class on the adaptation parameter. Furthermore, we 
investigate the effect of adaptation on the recruitment level, as well as on network topology. The 
theoretical predictions are compared with direct simulations of the full system. We identify two 
parameter regimes with qualitatively different bifurcation diagrams depending on whether nodes 
become susceptible frequently (multiple times in their lifetime) or rarely (much less than once per 
lifetime) . 
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I. INTRODUCTION 

Any society contains individuals who are carriers 
of an ideology or fad (e.g., a religious or political 
party afBliation) that they desire to spread to the 
rest of the society. Thus, for a given ideology, a 
society can be partitioned into a set of people that 
represent the ideology and want to spread it, and 
the complement of this set. For example, the ide- 
ology could correspond to the views of a particular 
political party, with the party members desiring to 
recruit new members to improve their positions in 
the government. Other areas of recruitment have 
been proposed as mechanisms for fads which appear 
as a rapid rise above some threshold, as in music 
[1] , management technologies [2] , economics [3] , and 
even science [4] . Slower recruitment based on social 
networks has been postulated in biology [5] and the 
spread of alcoholism [6] . Slowing the rise of a fad or 
eliminating the spread of an ideology has been also 
been proposed through the control of critical nodes 
in a social network [7] . 

In today's world where collisions between ideolo- 
gies can lead to radicalization of society, the prob- 
lem of existence and formation of terrorist networks 
or insurgency movements becomes important. Re- 
cently, mathematical modeling of various radical 
groups, such as terrorist organizations, has been 
done to explore their structure and dynamics [8- 
11]. In addition to dynamical approaches that mea- 
sure rates of attacks of radical groups [12], opera- 
tions research theories have also been applied to the 



formation of radical groups [13]. Another class of 
works [14-16] focuses on the recruitment dynamics 
of terrorists within a well-mixed population using 
compartmental models similar to those used for epi- 
demic spread. In [16], a systematic analysis of re- 
cruitment to Hezbollah has been done. Data from 
a political science discipline is used to generate a 
compartmental model of recruitment. Although the 
modeling is deterministic, it considers the various 
parameters which explicitly affect the success of re- 
cruitment to the radical's cause. Finally, a num- 
ber of recent publications discuss terrorist networks 
as optimal structures that optimize communication 
efficiency while balancing the secrecy of the net- 
works [17-19]. 

Not considered previously is the combination of 
the spread of radical ideas plus adaptive changes in 
social connections to improve recruitment success to 
the radical cause. Some work in this direction is 
presented in the papers studying voter models in 
which individuals make connections to influence oth- 
ers' opinions [20-22], as well as opinion dynamics 
models in which individuals are influenced by neigh- 
bors' opinions [23]. 

We develop a simple model of a society in which 
some of the members belong to a class that tries 
to spread an ideology. The ideology spreads as a 
result of contact of recruiting members with non- 
recruiting members. The recruiting members may 
improve their chances to recruit via network adapta- 
tion. The purpose of the adaptation proposed here 
is to improve the spread of the ideology, which is 
in contrast to network adaptation in epidemiologi- 



cal models [24, 25] where the purpose is avoidance 
of contact with the spreading members. In addition 
to adaptation, the connectivity within the society 
changes due to birth of new members and death of 
existing members. Using this model, we explore how 
the existence of stable recruiting classes depends on 
the adaptation and other parameters. We also con- 
sider the network topology of the recruiting class, as 
it may be important in accessing the quality of the 
communication channels in the resulting structure. 
Section II presents the model and a system of mean 
field equations describing its dynamics. Section III 
presents mean field analysis of the threshold for suc- 
cessful recruiting and its dependence on the network 
adaptation. These results are compared with simu- 
lations of the full system, and the network geometry 
of the recruiting class is considered. Section IV con- 
cludes. 



II. MODELING THE DYNAMICS OF 
RECRUITMENT 

A. Model 

Wc consider a social structure consisting of M in- 
dividuals represented by nodes in a network. An 
existing relation between any two individuals is rep- 
resented by a link between the two nodes. New in- 
dividuals join the society at a constant total rate fi, 
and the individuals in the network can leave the net- 
work via death at rate 6 per individual. When new 
individuals join the system, they arrive with a links, 
which are connected to a randomly selected nodes 
in the population. When individuals die, their links 
are removed. 
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FIG. 1. Schematic representation of node fluxes in and 
out of the network due to birth and death, described by 
the rates fi and 5 respectively. Fluxes between the three 
classes are described by the transition rates Ai, A2 and 
7. Link rewiring (not shown here) takes place at rate w. 

Some members of the society, recruiters (also re- 
ferred to as R- nodes), are carriers of an ideology 
that they try to spread to the individuals they come 
into contact with. The individuals that do not be- 



long to the recruiting class are divided into two 
groups: those who are susceptible to the recruit- 
ment (S-nodes) and those who are non-susceptible 
(N- nodes). We assume that individuals can sponta- 
neously change their state from non-susceptible to 
susceptible and vice versa (with rates Ai and A2 re- 
spectively). An individual in the recruiting class re- 
mains in that class until death [14]. 

A susceptible individual joins the recruiting class 
at a rate that is proportional (with proportionality 
constant 7) to the number of contacts it has with 
the recruiting class. In order to improve their re- 
cruiting capabilities, the recruiting individuals can 
rewire their links with rate w, abandoning a con- 
nection to a nonsusceptible individual in favor of a 
connection to a susceptible one. As the links are 
rewired, the network topology changes based on the 
current states of its nodes. A schematic representa- 
tion of the node dynamics is shown on Fig. 1. 

We perform Monte Carlo simulations of the above 
system following the continuous time algorithm de- 
scribed in [26]. We start with an Erdos-Renyi ran- 
dom network, which then evolves according to the 
rules of birth, death, and rewiring. The results pre- 
sented in the rest of the paper arc for the systems 
that have reached a steady state. 



Mean field 



To describe the dynamics of this system, we con- 
struct a mean field model for nodes and links as 
in [24]. The time evolution of the nodes of each 
type is described by the following rate equations: 



^tA/k = /i - AiWn + A2A/'s - (JA/k (la) 
dtMs = AiA/k - A2AA3 - 7AArs - 6J^s (lb) 
dtJVn = 7A/'rs - S^n, (Ic) 



Here the functions A/n = A/nC*), A/s = A/s(i), A/r = 
A/R(i) represent the number of nodes of each type. 
The process of recruitment is captured by the 7A/RS 
term, where the recruitment is shown to take place 
at a rate proportional to the number of links between 
the recruiting class and the susceptible class, Ars. 

In order to capture the rewiring process, we follow 
the evolution of the number of different types of links 
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present in the network: 
SjAnn = A2A/'sN + (TM 



A/n 



9t7VsN = o-^ 



A/k + A/'s + TVr 
2(Ai + 5)Mnn 
A/'s 



A/'n + A/'s + AAr 
(Ai + A2 + 2(5)AAsN + ZAiA/nn 



(2a) 

- 7A/'nsr + 2A2A/'ss 



-tA/ssr + AiWsN 
2(A2 + 6)Afss 



A/r 



StA/'ss 

a.A/'RN = 7AANSR + <.M_^^^_^^^_^^ 

- (Ai + 2^ + w)AArn + A2A/'rs 
i9tA/'Rs = -tA/rsr, + 7A/'ssR 

- (A2 + 2(5)A/'rs + (Ai + u;)AArn 
9tA/R,R, = 7A/'rsr, - S^A/'rr. 



(2b) 
(2c) 

(2d) 

(2e) 
(2f) 



where the terms A/xy correspond to the number of 
hnks connecting nodes from classes x and y, e.g., 
Ann corresponds to the number of NN links. The 
terms proportional to afi correspond to the influx 
of edges due to the birth of new nodes, where the 
probability for the new node to attach itself to a 
node from class X is proportional to the number of 
nodes in that class. The third order terms, Ansr, 
TVssr and A/'rsr, describe the formation of triples of 
nodes with an S-nodc at the center, and at least one 
of the edges terminating at an R-node. These terms 
describe the rate at which NS-, SS- and RS- links be- 
come NR-, SR- and RR- links respectively, due to the 
interaction of the central S-node with its neighbor- 
ing R-node. Note that our definition of RSR triples 
includes degenerate triples, i.e., RSR triples where 
both of the R-nodes correspond to a single R-node, 
by analogy with degenerate triangles. 

The resulting system of equations is not closed, 
as it contains higher order terms. Following earlier 
works in epidemiology [27], we introduce the closure 
based on the assumption of homogeneous distribu- 
tion of RS-links: 
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(3a) 
(3b) 

(3c) 



We thus obtain a system of nine mean field equa- 
tions, which we analyze. 



III. RESULTS 

A. Mean field recruiting threshold 

We first consider the bifurcation point of the re- 
cruitment model where the zero recruit (trivial) 
steady state becomes unstable, which we call the re- 
cruiting threshold. This is a transcritical bifurcation 
point, which is analogous to the epidemic onset in 
epidemic spreading models. It can be found analyt- 
ically as a function of parameters for the mean field 
model, and certain asymptotic limits have a simple 
form. 

We nondimensionalize Eqns. (1) and (2) to sim- 
plify the analysis. We introduce the dimension- 
less time variable t = St, where time is rescaled 
by the average node lifetime S~^. Further, the 
node variables are rescaled by the expected pop- 
ulation size at steady state, while the link 
variables are rescaled by the expected number 
of links at steady state, {fj,/S){a/2). Let x = 
[A^N, AAs , AAr, A/'NN2a-i , AAsn2(t-\ Mss2a-\ 
A/RN2cr~\7VRs2cr"\A/'RR2CT-i](5//~i denote the 9- 
dimensional vector of rescaled node and link state 
variables. Note that in the steady state the rescaled 
node variables sum to 1, and, therefore, in steady 
state they correspond to the probability for a node 
of a given type to exist in the system. Similarly, the 
rescaled link variables sum to 1, and correspond to 
the probability for an edge chosen at random to be 
of a given type, e.g., at steady state xs corresponds 
to the probability for a randomly chosen edge to be 
an RS link. 

We introduce the following rescaled parameters: 

(4a) 
(4b) 
(4c) 
(4d) 

We can now write the dimcnsionless equations of 
motion as 
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= Ai(5~i 


A2 


= A2(5^^ 
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± = F(x;A). 



(5) 



where A = [Ai, A2, VF, F, ct] is a vector of all the sys- 
tem parameters. (Recall that cr is a dimcnsionless 
integer.) The full system of dimcnsionless equations 
is given in Eqs. (AI). 

We can now find the trivial steady state so- 
lution, where the number of R nodes is zero. 
This restricts the state to be of the form Xq = 
[xo,i,xo,2,0,a;o,4,xo,5, 2:0,6,0,0,0], where the num- 
ber of links involving R nodes is zero as well. This 
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guarantees that Fi(xo, A) = for i = 3, 7, 8, 9. Since 
this subset of equations represents an invariant man- 
ifold, we concentrate on solving the equations for the 
rest of the 5 components, xq.i, ^0.2, 2;o,4, xo.5 and 
xo,6- 

The first observation is that the xi and X2 equa- 
tions are solvable when the number of recruiting 
nodes is zero, and they yield steady state values of 



xo,i = {A2 + l)D-\ 
xo,2 = AiL>-\ 



(6a) 
(6b) 



where Z) = Ai-|-A2-|-l. The nonzero link variables 
may be expressed in terms of xo,i, a;o,2, and they are 
given by the following: 



XOA = {A2 + lfD~\ 

xo,5 =2Al(A2 + l)^?-^ 



(be) 
(6d) 
(6e) 



Now that we have the full trivial solution, we can 
examine its stability. In order to do this, we lin- 
earize the vector field about Xq and examine where 
it has a one dimensional null space. That is, at the 
bifurcation point, there is only one real eigenvalue 
passing through zero. This is equivalent to examin- 
ing where the determinant of the Jacobian vanishes; 
i.e., we compute those parameters where 



det(2?xF(xo, A)) = 0. 



(7) 



The following relation describes the location of the 
bifurcation: 



W = -Ai 



2r 



-As 



[r(2-o-i)-l](Ai+A2 + l) 
1 2{T(7-^ + 1) 



r(2-cr-i)-i r(2-cr-i)-i 



(8) 



We next examine the limits of Eq. (8) when ei- 
ther the recruitment rate F or rewiring rate W is 
large. The minimum amount of recruitment required 
to maintain nonzero values of recruited population, 
when W approaches infinity, can be found by set- 
ting the least common denominator of Eq. (8) to 
zero, which yields a simple expression for F: 



F = 



1 



(9) 



If the recruitment rate is below this value, additional 
rewiring is not sufficient to enable spreading of the 
recruiting class. 



Similarly, we find the smallest amount of rewiring 
required for existence of a nontrivial steady state so- 
lution, conditioned on our ability to control F, with 
all the other parameters fixed: 



lim W 



-Ai 



1 

-1 



(2 



^)(Ai+A2 + l) 



For lower rewiring rates, the recruiting class cannot 
spread even if the recruiting rate is large. The value 
of this asymptote is greatest, meaning rewiring is 
most necessary, when Ai approaches zero. In this 
case few non-susceptible nodes become susceptible 
to recruiting, and the necessary rewiring value ap- 
proaches 
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(10) 



B. Comparison of mean field with simulations 

We next compare the mean field predictions for 
the spread of recruiters with simulations of the full 
stochastic network system. Thus, we compare the 
average size of the recruited portion of the popu- 
lation in statistical steady state to the solution of 
the mean field equations at steady state, which we 
solve exactly in the Appendix A. We assume that 
the parameters Ai and A2 (rates for gaining and 
losing susceptibility) and tr (determines average de- 
gree) depend on details of the society. On the other 
hand, we assume that the parameters F and W can 
be controlled by the recruiters, i.e., the recruiters 
may choose to be more or less aggressive in their re- 
cruitment, as well as in how quickly they rewire their 
links to susceptible members of society. Therefore, 
we investigate the recruiting effectiveness for a given 
choice of Ai, A2 and cr, while varying F and W. 

We distinguish two parameter regimes, corre- 
sponding to two qualitatively different behaviors of 
the system: small values of Ai (Ai <C 1) and large 
values of Ai (Ai ^ 1). The small Ai regime corre- 
sponds to the case where nodes become susceptible 
to recruiting only rarely, on average much less than 
once in their lifetime (where the average lifetime of 
a node is in the original time units and 1 in the 
dimensionless t units). The large Ai regime corre- 
sponds to the case where, on average, nodes become 
susceptible to recruitment at least once in their life- 
time. 
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FIG. 2. Direct numerical simulation of a stochastic net- 
work system with Ai ^ 1. In Fig. 2(a) we compare the 
result of the simulation (symbols) to the mean field so- 
lution given by Eq. (A18) (solid curves). Square (red 
online): strong rewiring, W « 40. Circle (green online): 
weak rewiring, W « 0.4. Fig. 2(b) shows the density plot 
representing the dependence of the fraction of recruited 
nodes in the population in statistical steady state on 
rescaled rates of rewiring, W, and recruiting, F. Vertical 
solid line (green online) corresponds to the minimum re- 
cruitment rate required to support a nontrivial solution, 
as given by Eq. (9). The solid curve (black online) corre- 
sponds to the recruiting threshold as predicted by mean 
field in Eq. (8). The simulations are performed on a sys- 
tem with = 10^ 5 = 1, o- = 10, Ai = 10, A2 = 100. 
The time averaged values are computed once the system 
reaches steady state regimes. (In all of the figures in this 
paper this corresponded to at least 105~^ units of time, 
in other words about 10 generations of nodes have died 
before we assume the system to be at the steady state.) 



In Figs. 2 and 3 we show the results of simulating 
the system in these two regimes. The density plots 
in Figs. 2(b) and 3(b) show the fraction of recruited 
nodes function of recruiting rate F and rewiring 
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FIG. 3. Direct numerical simulation of a stochastic net- 
work system with Ai ^ 1. In Fig. 3(a) we compare 
the result of the simulation to the mean field solution 
given by Eq. (A18). Square (red online): strong rewiring, 
W ~ 70. Circle (green online): weak rewiring, W « 0.2. 
Fig. 3(b) shows the density plot representing the depen- 
dence of the fraction of recruited nodes in the population 
in statistical steady state on rescaled rates of rewiring, 
W , and recruiting, F. Vertical solid line (green online) 
corresponds to the minimum recruitment rate required to 
support nontrivial solution, as given by Eq. (9) . Horizon- 
tal dashed line (red online) corresponds to the rewiring 
rate that guarantees existence of a nontrivial solution 
in the limit of large recruitment rate (Eq. 10). The 
solid curve corresponds to the recruiting threshold as 
predicted by mean field in Eq. (8). The simulations are 
performed on a system with jj. — lO'', 5 = 1, cr = 10, 
Ai = 0.01, A2 = 10. 



rate W in networks with Ai ^ 1 and Ai <g; 1 re- 
spectively. The black curves in the two figures rep- 
resent the location of the recruiting threshold as de- 
rived from the mean field equations and given by 
Eq. (8). Here mean field allows us to accurately pre- 
dict the onset of the stable nontrivial solution. In 



5 



Figs. 2(a) and 3(a) we compare the simulation re- 
sults to the mean field predictions. Even though 
there is some discrepancy between the direct simu- 
lations and the mean field near the bifurcation, the 
recruiting threshold and the asymptotic behavior for 
large F show excellent agreement with the simula- 
tions. 

We observe in simulations that the steady state 
trivial solution undergoes a forward transcritical bi- 
furcation in the recruiting rate F for all values of 
W for which a nontrivial solution exists. This is 
in contrast with epidemiological models where the 
purpose of the rewiring is avoidance of the nodes 
spreading infection [24, 25], which can undergo a 
backward transcritical bifurcation and exhibit bista- 
bility. Note that, unlike those models, the purpose of 
the rewiring in the recruitment model is attraction 
of the susceptible population by the recruiters. 

Wc would like to draw the reader's attention to an 
important difference in the system behavior in the 
two regimes. In the regime where Ai ^ 1, the non- 
trivial solution exists for all values of W as long as 
F is large enough. On the other hand, in the regime 
where Ai <^ 1, the nontrivial solution may fail to 
exist for any value of F unless rewiring is aggressive 
enough. In Fig. 3(b) there is a range of rewiring 
rates W for which only trivial solutions exist. The 
horizontal dashed line, given by Eq. (10), indicates 
the mean field level of rewiring that is sufficient for 
the nontrivial solution to exist for rapid recruiting 
(large F). Also, we can see that the rewiring level 
that is necessary for the emergence of the nontrivial 
solution can be very close to the predicted sufficient 
level of Eq. (10), for small values of Ai. 

Another important difference between the two 
regimes of Ai values is seen when we compare two 
systems with different values of Ai and A2 while the 
ratio Ai/(1 -I- A2) is kept fixed. Note that in the 
absence of recruitment, such systems have identical 
fractions of susceptible nodes. In other words, in 
the presence of recruitment, the pool of individu- 
als available for recruitment would appear to be the 
same. However, as we can see from Fig. 4, there is 
a significant difference in the size of the recruited 
population as well as the recruiting threshold. The 
difference appears to be caused by the difference in 
the dynamics in the two regimes. Thus, for the large 
values of Ai and A2 a node can become susceptible 
several times during its life, while as the value of 
Ai decreases, only some nodes will ever become sus- 
ceptible, with low likelihood of doing so more than 
once. 
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FIG. 4. Comparing systems witli different values of Ai, 
but same ratio of Ai/(1 + A2). Squares (red online), 
riglit axis: Ai — 10, A2 — 10^. Crosses (green online), 
left axis: Ai = 0.01, A2 = 9.001. Both of the simulations 
were performed with the following parameters: w = 10, 
At = 10^ (5 = 1, o- = 10. 
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FIG. 5. Direct numerical simulation. Density plot repre- 
senting the dependence of the fraction of the mean node 
degree in the R-subnetwork in steady state, 2A/'rr/A/r, 
on rescaled rates of rewiring, W , and recruiting, F. Solid 
line (green online) indicates the value of F where the de- 
gree is at a maximum for a given value of W , found using 
analytic expression in Eq. (B4). The white diamonds 
correspond to the location of the maximum in simula- 
tions for a given value of W . The simulations are per- 
formed with the following parameters: /i = 10^, (5 = 1, 
a = 10, Ai = 10, A2 = 100. 



C. Recruited subnetwork geometry 

We now investigate the structure of the portion of 
the network (later referred to as the R-subnetwork) 
consisting only of R-nodes and links between them. 
In particular, we are interested in the mean degree 
of its nodes when the system reaches a steady state. 
The mean field approach, in addition to predicting 
the fraction of recruited nodes (which corresponds 
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to the size of the R-subnetwork) , also provides in- 
formation about topology of the subnetwork in the 
form of mean degree of the nodes. The mean degree 
of a node in the subnetwork serves as a low order 
description of how well the nodes in the network are 
connected, which may be important for communi- 
cation within the established subnetwork. The den- 
sity plot in Fig. 5 shows the dependence of the mean 
degree on W and F. The nonmonotonic behavior 
as a function of F for fixed W is predicted by the 
mean field model. The solid line corresponds to the 
analytical prediction of the maximum's location as 
given by Eq. (B4) derived in Appendix B. Note that 
the analytic description of the maximum's location 
would allow the recruiting class to optimize connec- 
tivity within the subnetwork, if that happens to be 
an important goal for that class. The mean field ap- 
proximation fails to capture the nonmonotonic be- 
havior of the mean degree in the limit where Ai ^ 1, 
and we leave this issue to a future study. 

IV. CONCLUSIONS AND DISCUSSION 

We develop a toy model that describes recruit- 
ment of new members by an interested class within 
a society. The members of the recruiting class can 
improve their recruiting capabilities by via adapta- 
tion. Thus, they may choose to abandon their re- 
lations with those members of society that are not 
prone to recruitment. This model assumes that once 
a node joins the recruiting class, it itself becomes a 
recruiter and it remains a member of this class until 
death. In contrast to avoidance rewiring used to re- 
duce infection spreading in epidemic models [24, 25], 
network adaptation in our model promotes spread- 
ing. Additionally, the population is open with birth 
and death modeled explicitly, while most previous 
adaptive network models have been of closed popu- 
lations (e.g., [20-25, 28]). 

In this paper, we develop and analyze a mean field 
description of our model. Thus, we are able to ac- 
curately predict the onset of the stable nontrivial 
solution of the system at steady state. Furthermore, 
we are able to accurately predict the size of the re- 
cruited class for a given set of system parameters, as 
well as the mean degree of the subnetwork formed 
by the recruited members. We compare the predic- 
tions made by the mean field description with the 
direct simulations of our model. We generally ob- 
serve a good agreement between the model and its 
mean field approximation. 

Analyzing the mean field model, we find two pa- 
rameter regimes with very distinct qualitative be- 



havior. Thus, we show that in the society where 
the particular ideology is unpopular (perhaps corre- 
sponding to radical ideology) and individuals rarely 
become susceptible to it, adaptation is necessary in 
order to observe stable nonzero levels of the recruit- 
ing class. On the other hand, if the idea is suffi- 
ciently popular (e.g., recruitment into a moderate 
political party) , the adaptation improves the recruit- 
ing capabilities and may affect the ultimate topology 
of the recruited, but adaptation is not a necessary 
condition for the existence of a nonzero stable solu- 
tion. Furthermore, we speculate that if the model 
were changed to describe a society with two com- 
peting recruiting classes (think two party system), 
the adaptation may be a mechanism by which one 
competing class gets an edge over the other class. 

As evidenced by the results presented in Figs. 2(a) 
and 3(a), the mean field approximation has some in- 
accuracies in the parameter regime between the bi- 
furcation point and the asymptotic saturation. We 
predict that this inaccuracy is due to errors in the 
homogeneous closure assumption (Eqs. 3). In a sub- 
sequent work, we develop and analyze a new closure 
of the mean field equations that more accurately de- 
scribes the full system. 
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Appendix A: Exact solution 

In this appendix, we derive the exact solution to 
the system of equations in Eqs. (1) and (2) when the 
system is at the nontrivial steady state, i.e., when 
the left hand side of the equations is zero and the 
size of recruiting class is nonzero. We begin by de- 
riving the dimensionless equations, as defined in sec- 
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tion III A: 



respectively: 



dj:xi = 1 — Aia"! + A2X2 — xi (Ala) 

diX2 = Aixi - k2X2 - Fxs - X2 (Alb) 

d^x-i = Tzg - x-i (Ale) 

Xi 



dj:Xi = A2Xr, + 2 



dfX^ = 2- 



Xi+ X2 + X3 
X2 r^X5Xs 



- 2(Ai + l)x4 

(Aid) 



Xi + X2 + X3 X2 

(Ai + A2 + 2)2:5 + 2Aia;4 



2A2a;6 



(Ale) 



dixe = -2T^^ + Aixg - 2(A2 + 1)^6 (Alf) 



X2 



drX7 = r 



X5XS 



X3 



X2 X1+X2+ X3 

- (Ai + 2 + W)x7 + A2X8 
dixs = -r + 2<7-^xs^ + 2T 



yX2 J X2 

- (A2 + 2)xs + (Ai + W)xr 
d^xg = r + 2a^^xs^ - 2xg. 



(Alg) 

(Alh) 
(Ali) 



At steady state, the left hand side of the above equa- 
tions is zero. We proceed in our derivation by divid- 
ing all equations in the steady state by X2 and in- 
troducing a new variable Zi = Xi/x2-, obtaining the 
following system of equations: 



= 1/^2 - Aizi + A2 - zi 
= Aizi - A2 - - 1 
= Tzg - Z3 

= A22;5 + 2zi - 2(Ai + l)z4 
= 2- rzgzg + 2K2Z& - (Ai -I- A2 - 
+ 2A1Z4 

= -2TZ8ZQ + A1Z5 - 2(A2 + 1)Z6 
= rzszg + 2z3 - (Ai + 2 + W)zj - 
= -r (z| + 2(j-^Z8) + 2TZQZ8, 

-(A2+2)z8 + (Ai + VK)z7 
= r (z| + 2ct"^Z8) - 2z9. 



2)Z5 



A. 



(A2a) 
(A2b) 
(A2c) 
(A2d) 

(A2e) 
(A2f) 
(A2g) 

(A2h) 
(A2i) 



We have used the fact that in the steady state x'l + 
2^2+3^3 = 1, as can be shown by adding together 
Eqs. (la)-(lc). In the rest of the derivation we will 
find a closed equation for zg and express the other 
Zj 's in terms of zg. 

We can immediately express zi, Z3, and zg in 
terms of zg by solving Eqs. (A2b), (A2c) and (A2i 



zi = Ai'i(l + A2+rz8) (A3) 
Z3 = rzg (A4) 
Z9 = (r/2)(zg2 + a-izg). (A5) 

We substitute Eq. (A3) into Eq. (A2d), and Eq. (A4) 
into Eq. (A2g). The resulting two equations, to- 
gether with Eqs. (A2e), (A2f) and (A2h), form a 
closed system of five equations with five unknowns 
Z4-Z8. 

We solve Eq. (A2f) for Z5 in terms of zg and Z8: 



Z5 = 2Aii(l+A2+rz8)z6, 



(A6) 



and substitute the above result into Eq. (A2d) 
and (A2g), to solve for Z4 and Z7 in terms of zg and 
Z8 as follows: 

Z4 = [(Ai + l)Ai]-i[l + A2+rzg][l + A2Z6] (A7) 
Z7 = [rAr'(2rzg + 2(A2 + l))z6 + 2r + A2] X 
X (Ai +2 + W^)-izg. (A8) 

Substituting the expressions for Z4 and Z5 from 
Eq. (A6) and (A7) into Eq. (A2c), and solving for 
Zg we obtain 



Ai 



r(Ai + l)zg + (Ai+A2 + l)' 



(A9) 



Finally, substituting results of Eqs. (A9) and (AS) 
into Eq. (A2h) we obtain an equation for zg in a 
closed form: 



[(«ir)z| + (a2r 4- a3)zg + {ajV ^ + a5)]zg 
[r(Ai + l)zg + (Ai+A2 + l)] 



where 





(AlO) 
(All) 



ai = (Ai + l)(Ai + VK + 2) 

02 =2c7-i(Ai + l)(Ai+M^ + 2) 
-2(Ai +2)(VF + Ai) (A12) 

03 = 3(Ai + l)(Ai + VK + 2 + A2) + A2(VF + 1) 

(A13) 

04 = 2(Ai+A2 + l)(Ai + iy + 2 + A2) (A14) 
as = (Ai + A2 + l)[2(a-i - 2){W + Ai) 

+ 4(a-i-Ai)]. (A15) 

Note that the physically relevant solutions are posi- 
tive, and therefore finding a nontrivial solution of zg 
is a simple matter of solving a quadratic equation: 

) {aiT)zl + {a2T + a3)zg + [aiV-^ + 0,5) = 0. (AI6) 
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Solving Eq. (A2a) for X2, we can now express X2 
in terms of the newly found zg: 

X2 - Ai[(Ai + l)rz8 + (Ai + A2 + (A17) 

The rest of the original variables can be found using 
Xi — X2Zi. Thus, for example, X3 is 

X3 = X2Z3 = rAi28[(i + Ai)rz8 + (Ai + A2 + 1)]"^ 

(A18) 



Appendix B: Extremum of the mean degree in 
R-subnetwork 

In steady state, the mean degree of the nodes 
within the R-subnetwork is found by taking a ra- 
tio of twice the number of RR-links to the number 
of R-nodes in the subnetwork: 

^^^^^= A^ + '^r« + ^- 

where the values of A/rr and TVr are obtained by 
solving Eqs. (Ic) and (2f) in steady state. In this 
appendix we determine the value of F, Tm, that for 
a given rewiring rate will maximize the degree in the 



resulting R-subnetwork. We do this by maximizing 

Differentiating Eq. (AI6) with respect to F and 
evaluating the resulting equation at F,„, the value 
where extremum is attained, we obtain 



ai(-Z8)m + 02(2:8) 



rJai = 0. 



(B2) 



Note that the derivative of zg with respect to F 
evaluated at F^ is equal to zero because zg is an 
extremum there, and the are independent of F. 
Multiplying the above equation by F and subtract- 
ing it from Eq. (A16) evaluated at F„i allows us to 
solve for (zg),„: 



{Z8)r 



2a4 + Fas 
Fas 



(B3) 



Substituting the value of (zg)™ hito Eq. (B2) and 
solving for F,„ we obtain 

r,„ ~ [020304 - 2010405 

-I- (020304 -I- 01O5O4O3 — 02030504)^/^] (B4) 

/[05(0i05 - O3O2)], 

the location of the maximum for the given rewiring 
rate. 



[1] F. Meyer and A. Ultsch, in 
Advances In Data Analysis, Data Handling 
And Business Intelligence , German-Classification- 
Society (German-Classification-Society, 2010), pp. 
419-427. 

[2] J. Bendor, B. A. Huberman, and F. Wu, Journal Of 
Economic Behavior & Organization 72, 290 (2009). 

[3] M. A. Janssen and W. Jager, Journal Of Economic 
Psychology 22, 745 (2001). 

[4] E. Abrahamson, Scandinavian Journal Of Manage- 
ment 25, 235 (2009). 

[5] C. Grueter and F. Ratnieks, Animal Behaviour 81, 
949 (2011). 

[6] B. Benedict, SI AM News 40 (2007). 

[7] C. Kuhlman, V. Kumar, M. Marathe, 
S. Ravi, and D. Rosenkrantz, in 
Machine Learning and Knowledge Discovery 
in Databases. European Conference, 
ECML PKDD 2010. Proceedings (ECML, 2010). 

[8] A. Gutfraind, Terrorism as a mathematical problem , 
SIAM News (2009)^ 

[9] A. Gutfraind, Studies in Conflict and Terrorism 32, 
45 (2009). 

[10] N. F. Johnson, M. Spagat, J. A. Restrepo, 
O. Becerra, J. C. Bohrquez, N. Surez, 



E. M. Restrepo, and R. Zarama (2006), 
http:/ /arxiv.org/abs/physics/0605035. 

[11] A. Cherif, H. Yoshioka, W. Ni, and P. Bose (2009), 
arXiv/0910.5272. 

[12] N. Johnson, S. Carran, J. Botner, K. Fontaine, 
N. Laxague, P. Nuetzel, J. Turnley, and B. Tivnan, 
SCIENCE 333, 81 (2011), ISSN 0036-8075. 

[13] J. P. Caulkins, D. Grass, D. Feichtinger, and . Tra- 
gler, Comp. Oper. Res. 36 (2008). 

[14] F. Udwadia, G. Leitmann, and L. Lambertini, Dis- 
crete Dynamics in Nature and Society 2006, 85653 
(2006).' 

[15] K. Bentson, Master's thesis. Air Force Institute of 
Technology (2006). 

[16] L. B. Butler, Tech. Rep. 0MB No. 074-0188, 
SAMS Monograph, Fort Leavenworth, KS 66027- 
2134 (2011), contains many good references for re- 
cruitment, both stochastic and deterministic. 

[17] R. Lindelauf, P. Borm, and H. Hamers, Social Net- 
works 31, 126 (2009). 

[18] J. D. Farley, Studies in Conflict and Terrorism 26, 
399 (2003). 

[19] H. Bar-Isaac and M. Baccara, Working Pa- 
pers 06-07, New York University, Leonard N. 
Stern School of Business, Department of Eco- 



9 



nomics (2006), URL http://ideas.repec.0rg/p/ 
ste/nystbu/06-07.html. 
[20] I. J. Benczik, S. Z. Benczik, B. Schmittmann, and 
R. K. P. Zia, EPL (Europhysics Letters) 82, 48006 
(2008). 

[21] I. Benczik, S. Benczik, B. Schmittmann, and R. Zia, 
Physical Review E 79, 046104 (2009). 

[22] B. Schmittmann and A. Mukhopadhyay, Phys Rev 
E Stat Nonhn Soft Matter Phys 82, 066104 (2010). 

[23] P. Holme and M. E. J. Newman, Physical Review E 
74, 056108 (2006). 



[24] T. Gross, C. J. D. D'Lima, and B. Blasius, Phys 

Rev Lett 96, 208701 (2006). 
[25] L. Shaw and I. Schwartz, Physical Review E 77, 

066101 (2008). 
[26] D. Gillespie, Journal of Computational Physics 22, 

403 (1976). 

[27] M. J. Keeling, D. A. Rand, and A. J. Morris, Proc 

Biol Sci 264, 1149 (1997), URL http://dx.doi. 

org/10 . 1098/rspb . 1997 .0159. 
[28] T. Gross and B. Blasius, J R Soc Interface 5, 259 

(2008), URL http://dx.doi.org/10.1098/rsif. 

2007.1229. 



10 



